Helical Tubes in Crowded Environments 
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When placed in a crowded environment, a semi-flexible tube is forced to fold so as to make a more 
compact shape. One compact shape that often arises in nature is the tight helix, especially when 
the tube thickness is of comparable size to the tube length. In this paper we use an excluded volume 
effect to model the effects of crowding. This gives us a measure of compactness for configurations 
of the tube, which we use to look at structures of the semi-flexible tube that minimize the excluded 
volume. We focus most of our attention on the helix and which helical geometries are most compact. 
We found that helices of specific pitch to radius ratio 2.512 to be optimally compact. This is the 
same geometry that minimizes the global curvature of the curve defining the tube. We further 
investigate the effects of adding a bending energy or multiple tubes to begin to explore the more 
complete space of possible geometries a tube could form. 

PACS numbers: 81.16.Fg, 02.10.Kn,02.40.-k,05.20.-y 



I. INTRODUCTION 

The universe is a sparse and lonely place, with vast 
amounts of empty space. Most of us prefer to be in 
more dense regions, however, and it is precisely in these 
crowded environments that life occurs. These dense re- 
gions raise the question of how to pack things efficiently. 
Perhaps the most storied example is Kepler's conjecture 
that the face centered cubic (fee) lattice is the densest 
way to pack identical spheres. Though there was never 
much doubt that Kepler was right, a rigorous proof de- 
fied the centuries until Hales [l| developed a program of 
proof which established this long sought result. Free vol- 
ume theory [2J, for instance, uses the density of different 
lattices to establish their purely entropic free energy dif- 
ferences. As the fee lattice is the most densely packed, 
it allows for the greatest amount of translational entropy 
at any fixed volume fraction. When considering infinite 
volume systems, it is straightforward to pose optimality 
questions in terms of infinite periodic structures, whether 
it is the packing problem, the channel coding problem, 
the covering problem, etc. [3|. However, when consid- 
ering finite systems, or clusters, optimality can be more 
subtle. For instance, clusters of colloidal microspheres 
repeatably assemble into regular clusters that, in gen- 
eral, defy characterization through any packing principle 
[3, Q ■ Here we consider the optimal packing of a semi- 
flexible tube. Since the only degree of freedom for a sin- 
gle tube is its conformation, density is no longer a useful 
notion, as the tube's volume remains fixed. In this pa- 
per we consider the conformation of a solid, semi-flexible 
tube which maximizes the free volume for a hard sphere 
interacting only through its excluded volume. 

Other approaches have been proposed in defining the 
compact configurations of tubes. For instance, when 
studying knots made of solid tubes of fixed radius (in 
M 3 ), one can ask what geometry minimizes the length of 
the knot, i.e. what is the tightest knot? Though easily 
posed, this is a difficult question for which there is still 
no known minimizer for the simplest of knots, the trefoil 



Q . One approach is to minimize the global radius of cur- 
vature , which we will discuss in detail in the following 
section. Another approach is to minimize the convex hull 
of the tube's conformation. While both of these work for 
closed loops and knots, they are both problematic when 
considering finite, open tubes. In this case the straight 
tube is always optimal and so neither captures the notion 
of a compact conformation. 

Our approach, on the other hand, hinges on the 
Asakura-Oosawa depletion interaction [8|. We immerse 
our tube in a bath of hard spheres and maximize their 
entropy or, equivalently, their free volume. This gives 
conformations of the tube that capture the notion of com- 
pactness. The structure of this paper will be as follows. 
In section II we describe depletion volume theory, first for 
binary colloidal systems, and extend it to other systems. 
In section III we describe how to minimize the excluded 
volume of a free tube of finite thickness. We find that he- 
lices can arise for a range of tube lengths and support this 
with numerical integration of the free volume for spheres 
of varying radii. In the limit that the sphere radius is 
much smaller than the tube thickness, we are able to 
evaluate the free volume analytically and argue that this 
limit should generate those configurations found by min- 
imizing the global radius of curvature for closed knots. 
In section IV we consider the effect of additional interac- 
tions and study tube rigidity and tube-tube interactions. 
Section V summarizes our results. 



II. ENTROPY AND DEPLETION 
INTERACTIONS 

The depletion model of Asakura and Oosawa, first de- 
veloped to study polymer solutions, relies on the virial 
expansion for a mixture of large and small spheres 0]. 
Consider two large hard spheres of radius rj in a so- 
lution of smaller hard spheres of radius r s . The cen- 
ters of the two large spheres can get no closer than 
2ri apart while the small spheres can get no closer 
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FIG. 1: Excluded volume (green) for a straight and bent tube 
where the spheres cannot be. When the tube bends into the 
helical conformation there is overlap of the excluded volume 
between the layers of the helix and in the central core region 
of the helix. The excluded volume decreases by the overlap 
volume when the tube forms the helix. 



than ri + r s to the large spheres before they overlap. 
Thus we can enshroud each large sphere with a halo of 
thickness r s which excludes the small spheres. As the 
large spheres come closer, these excluded volume ha- 
los overlap and the free volume (Vf ree ) for the small 
spheres increases so that for N s small spheres, the free 
energy changes by AF = -k B TN s A\i\(Vf ree /V ) « 
-k B TN a /Wf ree /Vf ree = —k B TnV OV erlap where Vq is a 
characteristic volume for the small colloids, n — N s /Vf ree 
is the concentration of small spheres, and V over iap = 
AVf r ee is the overlap of the excluded volumes which is 
small compared to the free volume. This results in an ef- 
fective short range attraction between the large spheres 
at distances of order r s , an effect which has been seen 
in binary colloid systems Q and colloid-polymer systems 
[l(J. The free energy comes from purely geometric con- 
siderations: what is the arrangement of large spheres that 



maximizes the overlap of their halos and how large is that 
overlap volume? For two spheres the answer is obvious, 
the largest overlap occurs when the two spheres touch, 
and the overlap volume in this case is approximately 
litrir^ when r s -C r\. This result can be corroborated 
in detail through the virial expansion, for instance. 

The standard geometrical picture allows us to gener- 
alize, without developing the machinery of the virial ex- 
pansion, to other systems that interact via a hard-core 
interaction. For instance, a large sphere is attracted to 
a wall since its excluded volume overlaps with the wall's 
excluded volume. Thus the presence of boundaries com- 
plicates the problem of finding the configuration of large 
spheres which maximizes the overlap volume. Since the 
overlap of a sphere's and flat wall's excluded volume is 
roughly twice that for two spheres, approximately Aitrir 2 s , 
this attraction is large enough to drive two-dimensional 
crystallization along walls in place of three dimensional 
crystallization in bulk In this paper we replace the 
large spheres by a hard (but flexible) tube of radius t. In 
analogy to the halo around the large spheres there is a 
halo of excluded volume around the tube for all points 
within r s of the tube. In the previous example the large 
spheres had an effective attraction, while in our situa- 
tion the tube has an effective self attraction forcing it 
to bend. When the tube bends it is able to overlap the 
excluded volume from the two tube segments, leading to 
a decrease in the overall free energy. Thus we seek the 
shape of the tube which minimizes the free energy of the 
small spheres. We may also consider the optimal shapes 
of closed tubes and tubes tied into knots. We will discuss 
this in the following. 



III. FREE TUBES 

The space of possible configurations of the tube is dif- 
ficult to characterize, making it challenging to find the 
shape of the excluded volume (see, for instance, figure 
[T]). For an "ideal" tube (a tube that is perfectly hard, 
incompressible, but can bend flexibly) the configuration 



FIG. 2: If there is a section of tube where the excluded volume 
does not overlap with some other region of excluded volume 
then a deformation of the tube can increase the overlap. The 
move does not disturb the remainder of the chain, just as the 
center of a stack of coins can be shifted without moving the 
other coins. 
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space is additionally constrained by the tube thickness. 
Because the tube cannot bend and overlap itself, we 
can greatly reduce the configuration space by consider- 
ing shorter tube segments. At all segments of the tube 
the radius of curvature of the centerline must be larger 
than the tube radius, t, due to the tube thickness. As the 
tube gets longer, it becomes necessary to consider con- 
tacts from parts of the tube that are not "nearby" and 
then we use the global radius of curvature to account for 
this. 



A. Global Radius of Curvature and the Small 
Sphere Limit 

The radius of curvature at a point x on a curve C is 
the radius of the circle that passes through x and two 
adjacent points (in the continuum this can be defined 
through the obvious limiting procedure). Given three 



non-colinear points Xj 



1,2,3, we can construct the 



radius of the circle that passes through all three: 



r(xi,x 2 ,x 3 ) = 




(1) 



with 



w 



|x 2 - Xi| 

(x 3 - Xi) ■ (x 2 - Xi) 

|x 2 - Xi| 
|(x 3 - xi) x (x 2 - xi' 
|x 2 - Xi| 



(2) 



^2 f^ds , 



and where w ^ since the three points do not lie on 
a line. If we consider a continuous curve R(s), param- 
eterized by arclength s and choose Xi = R(s), x 2 
R(s + ds), and x 3 = R(s — ds), then u = ds + 
V = —ds/(l + -^Kds) and w = nds 2 / '(1 + ^nds), where k is 
the curvature of the curve. Then limd s ^o r ( x ij x 2,x 3 ) = 
1/k, recovering the radius of curvature. Using the 
fact that the area of the triangle is ^4(xi,x 2 ,x 3 ) = 
h\ (xj — Xj) x (xfe — Xj) | for and k distinct, we can 
rewrite dJ) as 



Xi -X 2 X 2 -X 3 X 3 -Xi 

r(xi,x 2 ,x 3 ) = — . 3) 

4A(xi,x 2 ,x 3 ) 

Gonzalez and Maddocks |7fl introduced the global ra- 
dius of curvature r g (x) to take into account contacts from 
points far away on the curve: 



r 9 (x) = min r(x, y, z) 



y.zGC 



(4) 



Just as the tube around C can be no thicker than the ra- 
dius of curvature, t < the global condition becomes 



t < A [CI = minr (x). 
xec 



(5) 
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FIG. 3: Allowed helix pitch and radius for a hard tube. He- 
lices with pitches and radii in the dashed region make unphys- 
ical configurations of the tube. When r a = t/10, only helices 
with pitches and radii in the shaded region have any overlap 
volume. 
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Helix Pitch (P/t) 

FIG. 4: Excluded volume as pitch increases above the min- 
imum allowed value P m in (dashed line). For small increases 
above P m in the excluded volume increases linearly with slope 
0.656 and plateaus as the turns separate enough so that there 
is no overlap of the excluded volume. We show our calculation 
for r s = t/10. 



In 0] it was shown that a closed curve of fixed length 
which maximizes the thickness A[C] has a constant value 
of r ff (x) = r* for all points x where it is curved and 
> r* when the curve is straight. We 



a value of r g (x) 

seek instead a curve which maximizes the overlap vol- 
ume V [C,r g /t] which depends on both the curve C and 
possibly the ratio of the radius of the deple ting spheres to 
the thickness of the tube. We have found [ll| that when 
r s /t — > 0, the tubes take on conformations for which the 
thickness of the tube t is precisely equal to the constant, 
global radius of curvature. Though this may not be sur- 
prising, it is certainly not obvious and does not hold in 
general - as the depleting spheres grow, the optimal tube 
configuration changes. 

We argue that as r s /t — > 0, if a smooth tube configu- 
ration exists which maximizes the overlap volume, then 
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r ff (x) = t for all points on the curve. To see this, con- 
sider a point x for which r g (x) ^ t: since the local radius 
of curvature cannot be smaller than t for a solid tube, 
r g cannot be smaller than t. Thus we consider the case 
that fg(x) = p > t. By assumption, the curve is smooth 
and so r g > p — e > t in a neighborhood of x. Thus 
there is an r s small enough so that there is no overlap 
of the excluded volume halo near x. Keeping r s fixed, 
we can move the curve in at x to achieve some overlap 
without reducing the overlap with some other part of the 
tube since there was no overlap at the outset. We de- 
pict this "stack of penny shift" in Figure [U Thus, we 
can increase the overlap volume if r g (x) > t and so that 
configuration cannot have the maximal overlap volume. 
Hence for small depleting spheres the ideal tube shape 
will have r g = t along the whole curve. We know no ar- 
gument to show that this is a unique configuration, and it 
is entirely possible that there are multiple configurations 
for which r g — t, only some (or one) of which maximizes 
the overlap volume. Note that this discussion required 
arbitrarily small depleting spheres. As the spheres grow, 
they begin to "see the forest for the trees" and, for in- 
stance, prefer straighter segments of tube to increase the 
local overlap volume - surfaces which bend away from 
the sphere have a smaller overlap volume than surfaces 
which bend towards the spheres. Thus, our approach al- 
lows us to both see the effect of finite depletor size and 
the optimal geometry of Gonzalez and Maddocks. 



B. Computation of the Excluded Volume 



An arbitrary point in the tube is at 

H(C, 0, s) = R(0 + a [cos 4, N(0 + sin B(£)] (8) 



for (j) £ [0,2vr], s G [0,t], and £ G [0,27m]. Not all val- 
ues of P and R are allowed [l2| when the tube is pre- 
vented from self-intersecting. The local curvature has 
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We start by restricting our study to helices as they 
are curves of constant curvature which can be chosen so 
that t = T g . Even when restricting our considerations to 
helices, there is a two parameter space of possible helix 
geometries. Aside from their length L, all helices can be 
described by their pitch P and radius R (see figure [TJ . 
The curve defining the centerline of the helix is 



R(£) = (i?cos(£),i?sin(£),P£/27r) 



(6) 



with £ G [0, 27m] where n is the number of helical turns 
the tube makes, n = , L In order to calculate 

^/P2 + (2x_R)2 

the overlap volume, we need to draw discs in the plane 
perpendicular to the curve's unit tangent T. A conve- 
nient basis for this plane is made by the normal N and 
binormal B unit vectors. With our parameterization, we 
have: 



T(0 

N(0 
B(0 



2tt 



[-i?sin£,i?cos£,P/27r] 



{2-kRY + P 2 

cos£, sin£,0] 
P 

(2nR) 2 + P 2 



[sin£,-cos£,27ri?/P] (7) 
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FIG. 5: Excluded volume as helix radius is changed for various 
sphere sizes. At small helix radius the geometry is constrained 
by the local curvature and the the layers of the helix do not 
touch. The big drop in excluded volume occurs when the 
helix turns start to have overlap between them. For small 
spheres r s — t/20 (A.) the minimum excluded volume occurs 
at the maximally compact helix with radius 0.86218. As the 
sphere sizes grow the minimum excluded volume occurs at 
larger radius, for spheres of radius r s = t/10 the minimum 
occurs at 0.88554 (B.) and at 1.01875 for spheres of radius 
r.=t/5(C). 
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an upper bound k — (p/ 2 J)1+r;i — j> an( i the distance 
between successive turns also must be greater than t. 
These constraints define the lower boundary of the ac- 
cessible region in Figure [3] Among all helices, the he- 
lix with P/R = 2.512 is the helix geometry where the 
thickness A[C] is equal to the local curvature and to the 
distance between successive turns, making it a compact 
helix [12j and the minimal excluded volume configuration 
when r s /t — ► 0. 

Fortunately, there is an upper bound to the relevant 
accessible region as well. If the helix is too rarefied, there 
will be no overlap of the excluded volume. Independent 
of tube configuration, an "ideal" tube with no self inter- 
sections will have volume V — irt 2 L. Similarly, if there 
is no overlap of the excluded volume, the total excluded 
volume will be K, xc i u dod = 7r(i + r s ) 2 L. We need not 
consider configurations for which the global curvature 
exceeds t + r s since an ideal tube of this larger radius 
can adopt these shapes. This consideration leads to the 
upper boundary of the region in Figure OH 

In the region between the two curves, we know of no 
simple formula for the total excluded volume and we have 
resorted to numerical integration. We take advantage of 
the screw symmetry of the helix to reduce the dimen- 
sionality of the integration. Recall that Pappus's theo- 
rem gives the volume of a surface of revolution as the 
area of the cross section multiplied by the length that its 
centroid traverses. Consider the straight line which lies 
at the center of the helix (the z-axis in our parameter- 
ization): the intersection of the helix with its excluded 
volume halo and any half plane beginning on this line has 
an area independent of the direction of the plane (i.e. the 
angle in the xy-plane) . The volume does not change if we 
slide each section up or down to be in registry with the 
cut through the irz-plane (for x > 0). Then, by Pappus, 
the volume is just the area of the cross section multi- 
plied by 27rr cen t ro id- To determine the area we pixelate 
the xz-plane and check if each point can be represented 
by H with s < t + r s . We can then calculate the area 
and the centroid. To properly compare the geometries 
we calculate the excluded volume per unit length, and 
find the optimal helix for a given sphere size. We find, 
not surprisingly, that for any helix radius, R, the mini- 
mum pitch always resulted in the largest overlap volume. 
This follows since increasing the pitch opens up the helix 
and allows gaps between successive turns of the tube. In 
Figure 2] we graph the overlap volume as we increase the 
pitch beyond its allowed minimum for fixed R and show 
how the excluded volume decreases as the pitch decreases 
at a given radius. For small increases in the pitch the ex- 
cluded volume increases linearly with P until successive 
turns no longer touch and then there is no longer an 
overlap region. The radius that minimizes the excluded 
volume depends on sphere size and we will denote it as 
R* in the following. 



C. Free Tube Results 

As we probe the space of possible helix parameters, 
the overlap volume goes through various regions of high 
and low overlap. This is depicted in Figures [5] A)-C) for 
sphere radii r s = t/20, t/10, and i/5, respectively. In 
all of these calculations, we set the pitch to its minimum 
allowed value for each R since this strictly minimizes the 
excluded volume. The largest amount of overlap occurs 
when layers of the helix lie one on top of the other. Thus, 
for small R when the layers can not achieve this geom- 
etry, the overlap is small. When the radius of the helix 
gets large enough, however, the local curvature does not 
prevent the turns of the helix from being near each other 
and the overlap volume greatly increases leading to the 
sharp decrease in excluded volume seen in the figures. 
The primary difference between different helices in this 
regime is the overlap in the central core region. This ac- 
counts for the rise in the excluded volume as R/t grows. 

We note that as the sphere size grows the optimal value 
of c = P/R decreases, though the general trend of ex- 
cluded volume versus R/t is unchanged. For instance, 
when r s = t/20 the minimum is almost exactly at the 
maximally compact helix with radius R* 05 = .86218i and 
P = 2.1658*, or c = P/R = 2.512. As the spheres get 
larger the optimal helix radius grows since the overlap 
of the excluded volume on the outside of the helix is in- 
creased by having a flatter contact. This effect competes 
with the overlap in the core and can dominate the en- 
tropy. When r s = t/10, we find fl^u = 0.885544* and 
P = 2.15444*, or c = 2.433, a smaller value than that 
for the most compact tube [ll|. As shown in Figure O c 
continues to drop as r s increases. 

As noted in earlier work by Maritan et al, [l7| . the 
optimal pitch to radius ratio is strikingly close to that 
for a-helices in proteins. This suggests that crowding it- 
self [HI can predict the geometry of proteins. Of course, 
hydrogen bonds stabilize the a- helix [14| , but our calcula- 
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FIG. 6: Ratio of the optimal helix pitch to helix radius (c = 
P/R) for various sphere sizes. As the spheres get smaller, the 
value of c approaches c* = 2.512 (dashed line). As the spheres 
get larger, c decreases because the helix radius is larger. 
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tion provides plausible reason for why that geometry was 
preferred in the first place. Our work is somehow dual 
to the use of site specific interactions which attempt to 
predict the secondary structure of proteins based on se- 
quence [l5|,[l6j|. It is in the spirit of Maritan et al. [l7l.ll8j 
who find this particular helical geometry as a recurring 
motif for a solid tube with a variety of interactions which 
serve to confine it and force a compact structure. Indeed, 
their work has shown the helix and other protein struc- 
tures as stable units which form purely on the basis of 
geometry. By constraining their study to tubes they were 
able to consider the intrinsic anisotropies of the system. 
In the protein, the anisotropy of the chain of amino acids 
comes from the side groups and the carbon backbone - 
both define the preferred direction The helix can 

even arise from tethered hard spheres if the bond length 
between neighboring spheres is different than the hard 
sphere radius [20j . 



D. Exact Analysis as r s /t — > 

We can gain further insight into our numerics by study- 
ing the small sphere limit where r s <C t. In this limit, 
overlap of excluded volume amounts to two segments of 
the tube touching. If the segments arc not parallel their 
contact will only be at a point and the overlap volume 
will depend on their relative orientation. If the tube seg- 
ments are aligned, however, then they will touch along a 
line. This is the case in the helix with turns lying one on 
top of another. Along a line of length £ the total overlap 
volume is the overlap area of two circles of radius r + ia 
distance 2t apart times the £. This gives an overlap 



V , - 1^.3/2^ 
v overiap — 3 



(9) 



We can also consider how this overlap volume changes 
when we pull the tubes apart to a distance d < r s : 



overlap 



4^2 



r 3/2 (1 _ a) 3/2 t l/2 £ 



., r s jH^(l-^a + 0(a 2 )) (10) 

where a = d/r s . In Figure 2] we have plotted the change 
in free volume per unit length as a function of the pitch 
P for fixed r s and R* s . This amounts to taking d — 
P — Pmin and we see that, indeed, the change in free vol- 
ume is linear in P and turns over when the tubes are no 
longer touching. The slope is m num « 0.66. We can use 
(fT0| to estimate the slope and find mtheory ~ Q.89£/L. 
What is the length I of the line of contact? By sym- 
metry, this line is itself a helix with the same pitch as 
the tube but with a different radius. Using the equation 
for the helix ([6]) we can consider where the piece of tube 
at £i = (xi = [R, 0, 0]) makes contact with another 
piece of tube centered at X2 = [i?cos(£), R sin(£), P£/2tt]. 
The contact line is halfway between these two centers, 



at x m id — ^[-R(l + cos(£)), i?sin(£), P£/27r] with a cor- 
responding radius p = i?|cos(£/2)| and helix length 
i = n^J P 2 + (2ttR cos(£/2)) 2 , where again n is the num- 
ber of turns in the helix. In general, the contact does not 
occur in a plane which contains the centerline of the helix 
- only in the case of stacked tori or, equivalently, infinite 
R, does this happen. Thus we find that i/L < 1. For 
fixed R and P we can find the value of £ which minimizes 
the distance |xi — X2I and find 



sin£ 



p 2 



(2irR) 2 



(11) 



Using our numerics for P/R when r s = t/10, we find 
roots at £ « ±1.77r so that £/L ss 0.91. This gives a final 
estimate of m t heory ~ 0.81, a reasonable, albeit sloppy, 
estimate for m„ um (« 0.66). Further correction is possi- 
ble by taking into account the reduction in overlap due 
to the local curvature. 

Armed with these analytic results, we can estimate at 
what length the helical conformation of the tube will give 
way to a larger superstructure (e.g. pleated sheets). For 
very short tubes, when L < irt, the tube can only form 
overlap of the excluded volume by bending into an arc of 
radius t. Once the tube is long enough to make contacts 
from distant points, it will from a nascent helix, provid- 
ing lines of contact instead of merely points of contact. 
The other possibility is that the tube bends back on it- 
self, making a pleat-like conformation. If we consider two 
helices next to each other, each of length L/2, we can de- 
termine whether the helix will bend back on itself by 
comparing the overlap volume of this configuration with 
that of a single helix. Splitting the helix into two seg- 
ments costs the contacts from one pair of turns, so the 
length of contact goes down by £\ ost > yj {2np) 2 + P 2 . 
However, the pair of helices now gains a contact length 
^gain < nP/2. Thus we expect the helix to be stable 
when nP/2 > y/ (2ivp) 2 + P 2 or, for the maximally com- 
pact helix with P/R w 2.512, we expect stable helices 
when n < 4.3. It is interesting to note that the average 
length of a-helices is 12 residues which corresponds to 
n = 3.3, well within our simple estimates [2lj . As the 
tube gets longer, the configuration that maximizes the 
excluded volume will surely get more complicated. In 
the limit of a very long tube, the problem becomes that 
of packing parallel tubes, an essentially two dimensional 
problem with the added cost of the hairpins necessary 
to bend the tube over. In this case the greatest overlap 
occurs when the tube bends back and forth to form an 
hexagonal lattice of cylinders similar to what is seen in 
hexagonally packed bundles of DNA. 



E. On the Knot but Not the Knot: Minimizing the 
Overlap Volume of a Trefoil 

We have argued that as r s /t — > 0, the depletion inter- 
action will drive the tube to its maximally compact con- 
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formation. Indeed, our result for the optimal helix reca- 
pitulates the geometry which minimizes the tube length 
for a given thickness. To further test our argument and 
make contact with the original motivation for developing 
the notion of global curvature, we will apply our analysis 
to the simplest knot, the trefoil. As with the helix, we 
choose a family of trefoil knots and minimize over the 
free parameters. There is a convenient class of trefoil 
knots, so called "torus knots", for which the centerline 
of the tube is embedded on a two-torus. Though this 
obviously does not span the whole space of such knots, 
it is convenient because of its analytic form: 

x(t) = [i? + &cos(67rt)]sin(47rf) 
y(t) = [i? + &cos(67rf)]cos(47rf) 
z(t) = &sin(67rf) (12) 

where R and b are the major and minor axis of the torus, 
respectively, and t G [0, 1] is the parameterization of 
the curve. Again, we define the tube by forming unit 
disks in the plane perpendicular to the curves tangent 
vector. Pieranski and Przybyl [22| found that among 
these embedded trefoils, the minimum global curvature 
occurs at R = 1.1158 and b = 1/2. As shown in Fig- 
ure [3 our numerical integration of the excluded volume 
shows the same trend. As R approaches 1.1158, the ex- 
cluded volume continuously decreases to its minimum at 
R = 1.1158. Moreover, R can be no smaller than 1.1158 
before the tube starts to self-overlap. Thus, as we argued, 
the depletion interaction drives the knot to its most com- 
pact structure as the depletors become small. 

Note that the ideal trefoil knot (i.e. the one that min- 
imizes the global curvature) cannot be tied on a torus. 
While there is no conformation proven to be the ideal 
trefoil knot, an upper bound on the minimum length of a 
trefoil knot has been found which cannot be achieved by 
a torus knot [13, [H| . This was found via the "shrink on 



no overlap" (SONO) algorithm, which is similar in spirit 
to maximizing the overlap volume. It is likely that these 
two approaches are equivalent as r s /t — > since our ar- 
guments related to global curvature should as well apply 
to the SONO algorithm. 

IV. BEYOND THE FREE TUBE: ADDITIONAL 
INTERACTIONS 

A. Tube Rigidity 

Real tubes, as opposed to the ideal tubes we have been 
considering, have a bending rigidity. The bending energy 
will compete with the entropy of the depleting spheres in 
selecting the ideal helix. As is standard with rigid rods, 
we introduce a bending energy which is controlled by t p , 
the persistence length (24|, the length scale over which 
the tube can bend with an energy of order ksT. The 
bending energy can be written locally in terms of the 
curvature of the tube, k(s): 

Feiastic = ^k B T£ p J K 2 (s)ds. (13) 

This can be specialized to the helix for which k is in- 
dependent of s and equal to k^ 1 = R [l + P 2 /(2ttR) 2 ] . 
Combining the bending energy with the entropy of the 
spheres allows us to calculate the change in free energy 
between a straight tube and a helix: 

AF « -k B TnV over i ap + kB ™P R2 (14) 

overlap 2 [R2 + p2 / ^ f] 2 K ' 

These two energies are minimized at opposite extremes. 
The bending energy is minimized for a straight tube with 
R/P = 0, while the excluded volume free energy is mini- 
mized by the compact helix. To investigate this balance, 
we find the densities and persistence lengths for which 
AF = between the straight tube and the compact he- 
lix for that sphere size. The competition is controlled by 
the relative size of n and £ p or, better, by the dimen- 
sionless quantities nr^ and £ p /t, related to the volume 
fraction of the spheres and the aspect ratio of the tube, 
respectively. We find it useful to define the control pa- 
rameter 9 = nrj? / (£ p /t). Again, for small r s /t, we can 

estimate V^ vor i ap /i ~ r Z J 2 t x l 2 and, for the ideal compact 
helix, k ~ l/t. Thus, we estimate the critical value of 
9 ~ (r s /i) 3 / 2 , which is confirmed numerically by Figure 
[HI Note that as r s Jt grows, the critical value of 9 drops 
below this value. This follows because as r s /t grows, the 
optimal helix is not as tightly wound and so the depletion 
effect need only overcome a smaller bending energy. 

We identify the control parameter as the ratio of two 
dimensionless parameters, <f> s — the volume 

fraction of the depleting spheres and £ p /t, the param- 
eter which controls the isotropic-to-nematic transition in 
the semiflexible tubes. This suggests that we have just 
probed one corner of a broader phase diagram where we 
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must also include the volume fraction of the tubes, 4>t- 
Roughly speaking, when <j> t <; 4t/£ p , the tubes should 
form a nematic phase. We sketch out a hypothetical 
phase diagram in Figure [5] It is an open and interest- 
ing question as to whether the tubes will form a nematic 
phase with equal numbers of left- and right-handed he- 
lices or, at a high enough concentration <p t , will form a 
cholesteric phase composed only of one type of helix. We 
can start to delve into that question by calculating the 
total excluded volume of two helices as they approach 
each other. We can get a sense of the difference in the 
overlap volume when two helices of opposite handedness 
are brought together versus two helices of the same hand- 
edness. 

In the case of two helices there is no special symme- 
try to exploit and simplify our calculation of the ex- 
cluded volume; we are now forced to do the full three 
dimensional numerical integration, a more computation- 
ally intensive endeavor. In this more complicated situ- 
ation there are many more degrees of freedom as well; 
the geometries of the respective helices, the distance be- 
tween the helices, the relative orientation of the helices, 
and the relative phase difference between them at the 
point of contact can all alter the excluded volume. For 
convenience we limit ourselves to the case where the two 
helices are both in the maximally compact configuration 
for isolated tubes and both have the same orientation so 
that overlap scales with length as opposed to contact at 
one point. This still allows us to explore the dependence 
of the excluded volume on chirality and distance between 
helices. All points in the first helix are defined as before 
by eq. [FJ and the points in a second helix of same hand- 
edness and no phase offset a distance D apart are defined 
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by 

H a (6 & a) 



R(£) + Dx + s [cos <j> N(£) + sin < 



■B(0] 
(15) 

for cp e [0, 2tt], s e [0, t], and £ e [0, 2im], and N, B, and 
R are defined by eq. [7] In this equation the chirality of 
the second helix is changed by replacing £ by — £, and a 
phase offset of if) radians is added to the second helix by 
replacing £ by £ + ip/(2ir). 

In Figure [TUR . we can see the overlap volume between 
helices of the same handedness and opposite handedness 
whose centers are at a distance 2R apart. The excluded 
volume is plotted as a function of the phase offset be- 
tween the two helices, and the minimum occurs at an 
offset of 7r when the two helices are just touching. As the 
offset either increases or decreases the helices stop touch- 
ing and the excluded volume increases. From this graph 
one can see that for all offsets, two helices of the same 
chirality have less excluded volume then for two of op- 
posite chirality. This result is misleading on the surface. 
Unlike the case of one tube, the local density of the tubes 
contributes to the entropy balance. In this case there is 
more overlap between two helices of the same handedness 
because they can not pack as densely and start to touch 
at further distances than two helices of opposite chirality. 
This is similar to the situation for hard spheres: the fee 
lattice is the densest packing of spheres so that at a given 
density of spheres, the fee lattice is the arrangement that 
gives spheres the most, not the least, free volume. Simi- 
larly helices of opposite handedness can move closer and 
will have more overlap volume at any given separation. 
One should also note that even ignoring chirality, the 
excluded volume dependence on offset shown in Figure 
[TUK appears strikingly similar to an anti-ferromagnetic 
XY model. This has been shown in dense packing of 
long chains to have complex crystallization which are ac- 
companied by spontaneous breaking of chiral symmetry 

In Figure HOB we show the excluded volume as enan- 
tiomeric helices are brought closer than 2R apart. Again, 
at every distance two helices of the same handedness have 
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we bring the helices closer. At any given distance, when the 
two helices have the same handedness there is less excluded 
volume, but the helices of opposite handedness are able to get 
closer, reducing the ultimate excluded volume. 



less excluded volume than two helices of opposite hand- 
edness, but eventually we get to the point where the two 
helices of equal handedness touch and can't move any 
closer. The two helices of opposite handedness can still 
move slightly closer and when they nearly touch the ex- 
cluded volume becomes lower than that for two helices 
of the same chirality at close packing. This shows that 
tubes of opposite handedness are favored as close-packed 
neighbors over two helices of the same chirality, but the 
difference is small. 

This brings us back to Figure [9] and our sketch of the 
phase diagram. The results here imply that for small <6* 
we will have an isotropic state. As <f>t grows beyond 4t/£ p , 
the tubes will align in a racemic state with roughly equal 



numbers of left- and right-handed helices. As <j>t grows, 
for sufficiently large <f> a , we expect a close-packed lattice 
of helices. Since it is impossible to decorate the close- 
packed triangular lattice in an anti-chiral fashion where 
every pair of nearest neighbors is of opposite handedness, 
we expect that the aligned, racemic nematic will give way 
to a cholesteric phase of helices of uniform handedness. 
In the middle of this phase diagram we can only speculate 
what would be seen - perhaps more complex double and 
triple helix bundles might form as the tube concentra- 
tion increases from the low concentration isotropic state. 
These complexes have been studied in the context of op- 
timal geometry and so we should expect them to appear 
as tubes are brought together [13] • From there the large 
number of degrees of freedom creates a vast array of pos- 
sibilities HI. 



V. CONCLUSION 

We have shown how the depletion interaction can drive 
helix formation and have estimated at what lengths the 
simple helix should give way to more complex tube con- 
formations. In addition to providing a local, force-based 
method for calculating the minimally compact curve, our 
approach also allows for direct comparison to experimen- 
tal systems with, for instance, wormlike micelles in a so- 
lution of spherical micelles. Our approach complements 
other purely geometrical approaches including the pio- 
neerin g w ork of Maritan, Micheletti, Trovato, and Ba- 
navar 17], the use of the global radius of curvature 0], 
and the use of overlap constraints The coincidence 
among the prior results and ours for the optimal helix 
strongly suggests that there is a primitive set of geomet- 
ric motifs that are the building blocks of long, tube-like 
molecules (29|. We have, in addition, found that maxi- 
mizing free volume finds the optimal trefoil knot and pro- 
motes the formation of paired, opposite-handed helices. 
Our model has the added advantage of another, indepen- 
dent parameter, the sphere radius r s . Thus the depleting 
spheres can act as surrogates for other interactions such 
as hydrophobicity and polymer-polymer interactions, or 
even as the side chains to amino acids when considering 
protein conformations. We have also considered the ef- 
fects of tube stiffness and compared these to the purely 
entropic depletion forces. Further work will more fully 
explore our proposed phase diagram (Figure [H]) and the 
higher order structures that are formed by the combina- 
tion of nematic interactions and depletion forces. 
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